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The  primary  objective  of  this  reseaurch  has  been  the  development  of  practical 
methods  for  the  compensation  and  control  of  multivariable  systems  in  order  to 
facilitate  the  design  of  controllers  for  complex  Air  Force  systems.  A  variety 
of  different  approaches  were  employed  to  accomplish  this  objective. 

In  particular,  the  authors  were  most  successful  in  developing  new  parameter 
adaptive  control  schemes  for  linear  multivariable  systems.  One  of  the  design 
algorithms  employs  an  adaptive  Luenberger  observer  which  v  (CONTINUED) 
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ITEM  #20,  CONTINUED:  automatically  adjusts  the  poles  of  a  given  system.  It 
is  felt  that  such  an  adaptive  controller  could  be  used  in  a  variety  of  aero¬ 
space  applications. 

A\ 

The  authors  trere  able  to  show  that  "multi-purpose"  controllers  can  be  designed 
«diich  simultaneously  perform  a  variety  of  control  functions.  In  particular, 
they  constructively  demonstrated  how'  to  build  a  controller  which  simultaneous! 
decouples,  places  poles  arbitrarily, ^rejects  disturbances,  insures  zero  error 
tracking,  and  is  robust  with  respect  to  peu'ameter  variations..  Additional 
studies  focused  on  the  implementation  of  such  controllers  using  digital  com¬ 
puters.  They  also  presented  a  new  and  straightforward  method  for  obtaining 
simple  low  order  models  of  systems  whose  dynamical  behavior  approximates  that 
of  more  complex,  higher  order  systems.  Such  low  order  models  can  be  used  in 
the  design  of  low  order  compensators  for  the  more  complex  systems. 

A  ccHnplete  new  resolution  was  also  presented  to  the  question  of  vdiat  changes 
occur  to  the  individual  transfer  matrix  elements  of  a  lineeu?  multivariable  ' 
system  under  constrained  output  feedback.  In  particular,  it  was  shown  what 
poles  become  controllable  and  observable  via  any  input/output  pair  tdien  con¬ 
straint  gain  output  feedback  is  applied  between  any  (i-th)  output  and  any 
(.1-th)  input.  The  changes  which  simultaneously  occur  to  the  numerator 
elements  of  the  transfer  matrix  were  then  determined  through  the  employment ‘R>f 
some  new  relationships  derived  from  an  approximate  relatively  right  prime 
factorization  of  the  system  transfer  matrix. 

Another  research  objective  was  the  study  of  parameterized  system  models  with  a 
view  towards  developing  compensation  and  control  techniques.  The  techniques 
considered  involved  the  methods  of  algebraic  geometry  and  revolved  around  thre« 
key  questions:  (i)  can  the  orbits  in  the  space  of  linear  systems  under  equiva> 
lence  via  the  action  of  an  algebraic  group  be  described  and  classified?, 

(ii)  what  spectral  structures  can  be  achieved  through  the  use  of  compensation?, 
and,  (iii)  what  are  the  essential  elements  required  in  extending  results  to 
domains  other  than  the  real  and  complex  numbers?  Question  1  was  resolved  for 
the  feedback  group  and  is  being  studied  for  systems  with  peurameters.  Results 
have  been  obtained  for  the  pole-assignment  problem  involving  parameters  using 
intersection  theory  and  some  preliminary  work  has  been  done  on  th  realiza¬ 
tion,  coprime  factorization  wd  trace  assignment  problems  for  systems  defined 
over  the  polynomial  ring  in  n-variables  over,  the  integers. 
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PRACTICAL  METHODS  FOR  THE  COMPENSATION  AND 
CONTROL  OF  MULTIVARIABLE  SYSTEMS 

1.  Research  Objectives 

The  primary  purpose  of  the  research  effort  has  been  and  continues  to  be 
the  development  of  practical  methods  for  the  analysis,^  compensation  and  control 
of  multivariable  systems  such  as  a  high  performance  aircraft,  helicopters, 
ballistic  missiles,  and  robots.  Success  in  this  development  facilitates  the 
design  of  complex  Air  Force  systems.  Four  critical  themes  form  the  foundation 
for  the  research;  namely:  (i)  the  question  of  parameter  variation  which  often 
involves  the  development  of  compensators  which  are  insensitive  to  that 
variation  (i.e.  which  are  robust),  (ii)  the  question  of  system  interaction 
which  leads  to  a  consideration  of  such  things  as  decoupling,  interconnected 
systems,  decentralized  control,  and  qualitative  properties;  (iii)  the  question 
of  stability  and  transient  response  under  compensation  which  typically  leads 
to  a  concern  with  pole-assignment  and  stabilization;  and,  (iv)  the  question  of 
computational  effectiveness  which  generates  a  search  for  constructive 
algorithms  with  a  view  towards  developing  automated  design  procedures. 

On  the  theoretical  level,  resolution  of  these  questions  has  required  no  e 
and  more  powerful  techniques  involving  algebra  and  algebraic  geometry. 
practical  level,  the  revolution  in  computer  hardware  has  created  an  excep¬ 
tional  opportunity  for  the  implementation  of  complex  designs,  especially, 
through  the  use  of  microprocessors.  Throughout  the  past  effort,  the  interplay 
between  the  theoretical  and  practical  levels  is  of  considerable  importance. 
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_2_  IIATIHW  J.  KERPER 

Chief,  Technical  Information  Division 

Section  two  of  the  report  contains  an  historical  survey  of  previous  work 
done  by  the  principal  investigators  under  Air  Force  sponsorship.  The  survey 
will  show  the  natural  evolution  of  the  research  in  the  broader  context  of 


system  science  and  engineering. 


2.  Historical  Survey 

System  science  and  engineering  are  concerned  with  the  modelling,  analysis 
and  control  of  phenomena  both  natural  and  man-made.  The  discipline,  whose 
formal  beginnings  go  back  at  least  to  Watt  and  Maxwell,  derives  much  of  its 
motivation  from  control  engineering  problems.  Before  World  War  II,  system 
design  and  analysis  were  primarily  an  art  and  the  key  concern  was  the 
regulation  of  simple  scalar  systems.  During  and  after  the  war,  techniques 
based  on  complex  variable  theory  and  the  Laplace  transform  were  developed  and 
applied  to  single-input,  single-output  systems  represented  by  a  rational 
function  called  the  transfer  function.  The  fire  control  problem  was  a 
particularly  important  application  in  this  era.  The  theory  of  servomechanisms 
developed  rapidly  from  the  end  of  the  war  through  the  Fifties  and  time— domain 
techniques  began  to  be  applied.  Most  controllers  and  sensors-  were  analog 
devices.  The  beginning  of  the  space  age  with  its  guidance  and  control  requir¬ 
ements  and  the  representation  of  transfer  functions  via  linear  constant 
coefficient  differential  equations  led  to  an  interest  in  state  space  methods. 
The  introduction  of  the  Kalman  filter  and  LOG  (linear-quadratic-gaussian) 
design  methods  had  considerable  impact  on  the  treatment  of  multivariable 
systems.  Frequency  domain  methods  using  polynomial  matrices  were  also 
developed.  However,  the  increasing  complexity  and  scale  of  the  engineering 
and  economic  problems  considered  required  greater  sophistication  on  both  the 
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theoretlcal  and  practical  levels.  Resolution  of  the  theoretical  issues  has 
required  more  and  more  powerful  techniques  involving  algebra  and  albegraic 
geometry.  On  the  practical  level,  implementation  of  the  more  complex  design 
procedures  has  led  to  a  greater  concern  with  the  use  of  digital  computation. 

The  broad  evolution  of  the  field  as  well  as  the  four  crucial  themes  have 
impacted  the  work  of  the  principal  investigators.  Early  work  on  decoupling 
and  system  structure  (1),  (2)  paved  the  way  for  considerable  later  efforts  on 
non-interaction  and  decentralized  control.  In  addition,  results  were  aimed  at 
providing  a  practical  design  approach.  To  illustrate,  the  necessary  and 
sufficient  condition  of  (1)  led  to  an  easily  implemented  numerically  stable 
algorithm  for  pole  placement  which  was  used  in  the  manual  control  of  a 
hovering  helicopter  (13).  Considerable  research  was  devoted  to  the  character¬ 
ization  of  those  "parts"  of  a  system  which  remain  invariant  and  those  "parts" 
which  can  be  freely  adjusted  in  various  ways.  Initial  general  results  occur 
in  (4)  and  the  principal  investigators  obtained  results  giving  considerable 
insight  regarding  the  effects  of  dynamic  compensation  (5)  and  feedback(6).  In 
(5),  the  "interactor"  of  a  transfer  matrix  was  introduced  and  represents  a 
quantitative  as  well  as  qualitative  measure  of  system  coupling.  This  connec¬ 
tion  was  exploited  in  (7)  and  has  been  related  to  invertibility.  In  view  of 
the  complexity  of  such  design  methods  as  dynamic  compensation  and  in  view  of 
the  desirability  of  treating  systems  with  parameter  variation,  research  on 
adaptive  control  ((8),  (9),  for  example)  and  on  pole  assignment  with  para¬ 
meters  (10)  was  undertaken.  While  initial  work  on  adaptive  control  involved 
the  model  reference  approach,  a  new  type  of  approach  was  developed  which 
resulted  in  a  combination  of  on-line  identification  and  control  (pole  place¬ 
ment).  Further,  the  analysis  of  adaptive  control  resulted  in  an  interest  in 
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low  order  models  as  well  as  model  matching,  tracking,  and  regulation.  A  new 
and  straightforward  method  for  obtaining  simple  low  order  models  of  systems 
whose  dynamical  behavior  approximates  that  of  more  complex  systems  was 
developed  (11).  It  was  also  shown  that  "multipurpose"  controllers  can  be 
designed  which  simultaneously  perform  a  variety  of  control  functions.  In 
particular,  it  was  shown  that  a  controller  could  be  built  which  simultaneously 
decouples,  places  poles,  rejects  disburbnces,  insures  zero  error  tracking  and 
is  robust  with  respect  to  parameter  variations  (12).  As  the  research 
progressed,  it  became  clear  that,  on  the  theoretical  level,  the  powerful 
techniques  of  algebraic  geometry  had  a  crucial  role  to  play  and  a  broad 
summary  was  given  in  (13).  It  also  became  clear  that,  on  the  practical  level, 
the  most  important  issue  is  the  question  of  actual  implementation  of  the  final 
design  and  that  digital  computation  will  play  a  key  role  in  that  implementa¬ 
tion. 

Summarizing  then,  significant  progress  has  been  made  with  respect  to  the 
development  of  controllers  for  multivariable  systems,  even  when  parameters 
vary  or  are  unknown.  Numerous  techniques  are  now  available  for  achieving  a 
variety  of  design  goals  and,  in  particular,  classical  scalar  techniques  such 
as  root  locus  and  the  Nyquist  stability  criterion  have  been  extended  to  the 
multivariable  case  (14)  -  (19).  Note  that  (18)  and  (19)  use  algebraic 
geometry  in  an  essential  way.  The  robustness  properties  of  LQG  designs  have 
been  demonstrated  (20),  (21)  and  techniques  have  been  developed  for  achieving 
a  diversity  of  design  goals  simultaneously  using  a  single  controller  (albeit 
of  rather  high  order)  (12),  (22),  (23).  Adaptive  control  has  also  received 
new  impetus  in  light  of  recent  theoretical  developments  (24)  -  (26)  coupled 
with  a  differential  operator  approach  (27)  -  (29).  Hoiraver,  unlike 
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classical  control  and  compensation  techniques  which  often  require  little  more 
than  pencil  and  paper  computations,  many  of  the  modern  multivariable 
procedures  require  extensive  computation.  Indeed,  computers  and  computer 
aided  design  algorithms  are  frequently  essential  in  achieving  a  final  design. 
3.  Applications  of  Algebraic  Methods  in  Systems  Science  and  Engineering 

As  noted  in  the  introduction  and  historical  review,  there  have  been  four 
critical  themes  in  the  research  effort;  namely:  (i)  the  question  of  parameter 
variation  and  the  devleopment  of  compensators  which  are  insensitive  to  that 
variation  (i.e.  which  are  robust);  (ii)  the  question  of  system  structure  and 
interaction  which  lead  to  a  consideration  of  such  things  as  interconnected 
systems,  decentralized  control,  and  qualitative  properties;  (iii)  the  question 
of  stability  and  transient  response  under  compensation  which  results  in  a 
consideration  of  pole  assignment;  and,  (iv)  the  question  of  computational 
effectiveness  which  involves  a  search  for  constructive  algorithms  with  a  view 
towards  developing  automated  design  procedures.  On  the  theoretical  level, 
resolution  of  these  questions  has  involved  the  use  of  the  techniques  of  modern 
algebra  and  algebraic  geometry. 

The  particular  conceptual  context  of  the  work  revolved  around  three  key 
questions,  namely: 

Question  1  -  Can  the  orbits  in  the  space  of  linear  systems  under  equivalence 
via  the  action  of  an  algebraic  group  be  described  and  classified? 

Question  2  ■*  What  spectral  structures  can  be  achieved  through  the  use  of 
compensation? 

Question  3  -  What  are  the  essential  elements  required  in  extending  results  Co 
domains  .other-  than  the  real  and  complex  numbers? 
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A  svumary  of  various  results  and  problems  is  given  in  (13).  In  terns  of 
the  critical  themes,  Question  1  relates  to  the  issues  of  parameter  variation, 

system  strucure,  and  stability;  Question  2  relates  to  these  issues  and  the 
issue  of  computational  effectiveness;  and  Question  3  relates  particularly  to 
the  issue  of  system  structure  and  more  importantly,  to  the  problem  of  combin¬ 
ing  digital  and  continuous  systems. 

Some  of  the  significant  results  obtained  during  the  grant  period  are  des¬ 
cribed  and  summarized  in  the  sequel. 

1.  The  Inter actor  and  Decoupling 

Given  any  proper  rational  transfer  matrix  T(s),  a  special  lower  left 
triangular  matrix,  C^Cs),  called  the  interactor  has  been  defined  and  shown  to 
be  (together  with  the  rank  of  T(s))  a  complete  invariant  under  dynamic  compen¬ 
sation  (5).  Results  on  decoupling  and  pole  placement  were  obtained  (7).  The 
main  theorems  are: 

Theorem:  A  system  characterized  by  a  nonsingular,  proper,  rational  mxm  trans¬ 
fer  matrix  T(s)  can  be  decoupled  via  state  feedback  if  and  only  if  the  inter¬ 
actor  Ct(s)  is  diagonal. 

Theorem:  Let  T(s)  be  a  proper  rational  pxm  transfer  matrix  of  full  rank  p. 

Let  P(3)  be  a  proper  rational  pxp  diagonal  transfer  matrix  such  that 

is  proper.  Then  there  is  a  dynamic  compensator  Tq(s)  such  that 
T(8)Tc(s)  -  D(s). 

Theorem:  A  system  characterized  by  a  nonsingular,  proper,  rational  mxm 
transfer  matrix  T(s)  can  always  be  triangularly  decoupled  with  all  closed  loo] 
poles  assigned  using  state  feedback. 

These  results  relate  to  Question  2. 


2.  Invariants  and  Canonical  Forms  under  Feedback 

A  complete  resolution  of  the  problem  of  determining  invariants  and 
canonical  forms  under  feedback  for  linear  systems  characterized  by  proper 
rational  transfer  matriices  was  obtained  (6).  The  invariants  were  determined 
in  the  frequency  domain  and  consist  of  the  Kronecker  set  of  controllability 
indices  together  with  a  canonical  form  for  the  numerator  of  the  transfer 
matrix  under  the  action  of  a  stabilizer  subgroup  of  the  unimodular  group  of 
polynomial  matrices. 

Let  T(s)  be  a  proper  (full  rank)  pxm  transfer  matrix  so  that  T(s)  “ 
Rt(s)P<j>~1(s)  where  R<r»  Pt  are  relatively  right  prime  polynomial  matrices  with 
Pf  column  proper.  The  set  »  {3l(Px)»  •  •  •  cb^PT^  )  called  the 

Kronecker.  set  of  T.  If  8l  21  ^  1.  •  •  >  3  m>  then  Pj  is  properly 

indexed.  Let  n  »  n^  “  deg  det  and  let  g  be  an  element  of  GL(m),  i  be  an 
element  of  (  =  nxm  matrices).  Call  (f,g)  a  state  feedbck  pair. 

Set 

PTf,g  =  8"^  {Pt  -  »  RTfg“  Rt.  Tf,g=  Rxj^g  PXf^g 

Then  Tj  is  equivalent  to  T  under  state  feedback  if  there  are  state  feedback 
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pairs  <f,g),  with  Tf^g  »  Ti,  Tif^g  “  T.  This 

notion  defines  an  equivalence  relation.  If  0(T)  is  the  orbit  of  T  (i.e. 
equivalence  class),  then  there  is  a  Tj  in  0(T)  with  properly  indexed. 

Thus,  in  determining  invariants  and  canonical  forms  it  is  enough  to  deal  with 
properly  indexed  systems.  So,  let  W  ■{  p|  P  properly  indexed  with  Kronecker 
set  3}  and  let  =  ^^*^3  )  “  |  W^U  ■  ,  U  unimodular)  be  the  stabilizer 

of  .  Let  Xg  *  {  r(3^(R)  _<  £}.  Then  is  stable  under  the  action  of  S 

on  the  right.  The  main  results  are: 

Theorem  A  canonical  form  for  the  action  of  on  exists  (and  is  denoted 

^  -■  -  ■  —  ■  ■  ,  —  —  ■  -  d  '  0  <1.  ■  ■■  —  — 

by  R^). 

Theorem  A  complete  system  of  invariants  for  equivalence  under  state  feedback 
is  given  by  (Re , 9  ) . 

The  proof  of  the  first  theorem  relys  heavily  on  algebraic  methods. 

3.  Pole  Assignment  with  Parameters 

Consider  the  system  x(p)  ■>  (A(p),  B(p))  where  A,B  depend  algebraically  on 
the  parameter  p=  (pj  ,  ...  p^).  Let  Pn^s)  x(s)  :  x  a  monic  polynomial 
of  degree  n)  and  let  ^(p):  k®”  “*■  P„(s)  be  the  affine  map  defined  by 

”x(p)  -  det  I  si  -  (A(p)  +  B(p)F)  1 

-  a  (A(p)  +  B{p)F*  Ojj(s) 

where  a(A(p)  ♦  B(p)F)  -  (la  i.(A(p)  +  B(p)F)...)  is  the  vector  of  character¬ 
istic  coefficients  and  •••  s  D*  The  problem  is  to 

determine  conditions  which  ensure  that  there  is  a  subspace  ^  of  k®^  such 
that 


^x(p)^^)l^  “  Pj(8)'Ks,p) 

for  some  <1/  (a, p)  of  degree  n  -  £  in  st  In  such  a  case,  it  isay  be  said  that 
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poles  can  be  assigned  independently  of  p  .  Results  on  this  problem  were 
derived  in  (10)  for  the  case  A(p)  =  Aq  +2  p£A£,B(p)  =  Bq.  Recently,  more 
general  results  have  been  obtained  using  the  methods  of  algebraic  geometry* 

In  particular,  there  are  the  following  theorems: 

Theorem  Let  A(p)  =  Aq  +  AjCp)  +  ...  +  Aj.p)  where  A£  is  homogeneous  of 
degree  i  and  let  B(p)  =  Bq.  Suppose  that  A(p)  is  non-degenerate  with 
respect  to  p  Then  n  poles  can  be  assigned  independently  of  p  if  and 
only  if  A£(p)  =  0  for  i  ^  1  and  (Aq,Bo)  is  controllable. 

(This  theorem  shows  that  the  problem  has  content  in  the  sense  that  if  there  is 
true  dependence  on  the  parameter,  then  completely  independent  pole  assignment 
is  not  possible.) 

Theorem  Let  A(p)  =*  Ap  Ai(p)  +  ...  *  A  (p)  where  A£(p)  is  homogeneous  of 
degree  i  and  let  B(p)  =  B^.  Then  11  x(p)^^^  “  ^x(p)  ®  subvariety  of  the 

affine  space  of  polynomials  Pn,vn^®>P^  (monic  in  s)  and  so  is 

(p)(F^)  for  any  subspace  of  k™n.  Moreover,  there  is  an 

Fj  such  that  Vx(p)(F£  )  =  P£(s)i|>  (s,p)  if  (i)  there  is  an  fg  in  F  such 

- - 1 - -  i 

that  ^(p)  (fp)  =  3  (s,p)  and  (ii)  the  varieties  Vx(p)^^ie,^ 

Pjjfsl'Ksjp)  in  Pn,vn^®»P^  transversal  at  ^x(p)^^o^*  (This  theorem 

provides  a  necessary  condition  for  the  problem.) 

Calculation  of  the  maximum  I  for  which  a  solution  exists  were  derived  in 
(10)  for  the  case  A(p)  =  Ag  +  1^P£A£,  B(p)  =  Bg. 

A.  System  Structure 

Consider  the  set  of  linear  system  (A,B,C)  with  A  nxn,  B  nxm  and  Cpxn. 
Then  GL(n)  acts  on  (A,B,C)  via 


♦This  means  that  R(Bgl ^  Ker  (A£(p)  )  ■  (0)  i  >  0. 


1 
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g*  (A,B,C)  =  (gAg”^-,gB,Cg~^,  gJGL(n) 

This  action  defines  state  space  equivalence  and  it  is  natural  to  ask  whether 
an  algebraic  quotient  exists.  In  (4),  it  was  shown  that  the  field  of  invar¬ 
iants  =  (k(A,B,C))®^^*^^  is  the  same  as  the  function  field  K(CB,  ... 

Ca2"~^B)  of  the  Hankel  matrices.  If  ^  (A,B,C)  :  (A,B)  controllable, 

(A,C)  observable  then  the  following  theorem  was  proved  in  (13). 

n 

Theorem  A  geometric  quotient  state  space  equivalence  exists. 

Moreover,  p^n,  is  a  rational,  smooth,  irreducible  quasi-projective  variety  of 
dimension  nCm+p). 


This  theorem  has  applications  to  system  identification  and  to  the  determin¬ 
ation  of  generic  properties. 

Increasing  interest  in  image  procesing  and  in  the  interconnection  of 
digital  and  continuous  systems  and  in  systems  with  fixed  delays  has  led  to  a 
consideration  of  systems  involving  transfer  matrices  whose  entries  are 
rational  functions  in  n-variables.  A  key  result  that  was  obtained  is  the 
following  lemma: 

Lemma  Let  R  be  a  Noetherian  integral  domain  of  characteristic  0  and  let  Z 
be  the  integers.  Then  the  realization  and  coprime  factorization  problems  are 
solvable  over  R[xi,...,X(,  ]»  Z  [xi,..,Xn  ]  R  if  and  only  if  they  are 
solvable  over  Z  [xi,...,x„  ]  and  R  is  a  regular  extension  of  Z. 


This  lemma  exposes  the  essential  algebraic  nature  of  certain  basic  system 
problems  and  has  considerable  implications  for  system  structure. 
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5.  Adaptive  Control  and  Model  Reduction 

The  main  thrust  of  our  work  in  this  area  has  been  the  development  of 
practical  algorithms  which  can  be  used  to  adaptively  control  multivariable 
systems.  At  the  outset  of  this  work,  considerable  progress  had  been  made  in 
designing  globally  stable  adaptive  controllers  for  unknown  continuous  time 
scalar  systems  (24)  -  (26).  Little  progress  had  been  made  in  extending  such 
designs  to  the  multivariable  case.  This  was  remedied,  in  part,  by  the 
multivariable  adaptive  control  scheme  detailed  in  (28).  In  particular,  in 
(28)  an  equation  is  derived  chat  can  be  used  to  identify  a  set  of  controller 
parameters  which  completely  assign  all  poles  and  zeros  of  a  linear  multivari¬ 
able  systems  represented  by  a  (pxm)  strictly  proper  transfer  matrix  T(s)  when 
p  ^  m  and  only  input-output  data  is  available  for  measurement.  This  equation 
is  used  Co  derive  a  parameter  adaptive  control  scheme  for  linear  multivariable 
systems . 

The  control  structure  contains  an  adaptive  Luenberger  observer  which 
assigns  as  poles  of  Che  closed  loop  system  Che  zeros  of  the  unknown  system 
(plant),  and  possibly  some  additional  poles  specified  by  the  designer.  The 
new  overall  system  zeros  and  the  remaining  poles  are  obtained  by  use  of  a 
fixed  precompensator  which  can  be  arbitrarily  specified  by  Che  designer. 

Since  Che  incorporated  parameter  identifier  directly  estimates  the  control 
parameters  without  explicitly  identifying  a  parameterized  model  of  Che  unknown 
plant,  the  structure  can  be  classified  as  a  direct  adaptive  controller. 

Since  Che  plant  zeros  are  tissigned  as  closed  loop  poles,  to  construct  Che 
controller  one  snisc  be  assured  Chat  Che  plant's  open  loop  zeros  lie  strictly 
in  Che  left  half  plane.  In  addition,  to  complete  Che  design  one  must  know,  a 
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priori)  an  upperbound  on  the  plant's  observability  index,  as  well  as  know  the 
structure  of  the  interactor  matrix  (5)  associated  with  T(s).  In  many  cases, 
this  latter  requirement  is  essentially  equivalent  to  knowledge  of  the  smallest 
relative  degree  in  each  row  of  T(s).  To  assure  proper  performance  additional 
information  about  the  plant  transfer  matrix  structure  is  also  needed.  This 
information  is  analogous  to  information  on  high  frequency  gain  necessary  for 
adaptive  control  of  scalar  systems  (24).  It  might  be  noted  that  the  design 
presented  in  (28)  is  more  general  than  those  presented  in  earlier  reports 
since  it  is  applicable  to  a  larger  class  of  multivariable  systems. 

More  recently,  the  restrictions  associated  with  this  direct  adaptive 
controller  have  been  partially  overcome  by  the  development  of  a  new  adaptive 
control  structure  for  MIMO  systems.  It  is  applicable  to  the  class  of  systems 
characterized  by  proper  full  rank  transfer  matrices,  and  can  be  implemented 
with  only  a  priori  knowledge  of  the  system  controllability  indices.  The  adap¬ 
tive  scheme  is  based  upon  a  fixed  observer  structure,  and  arbitrarily  assigns 
the  closed  loop  system  poles.  The  algorithm  can  be  viewed  as  a  generalization 
of  the  SISO  scheme  for  direct  adaptive  pole  placement  discussed  in  (30).  A 
preliminary  version  of  this  technique  was  reported  at  the  Twentieth  Decision 
and  Control  Control  Conference  (31).  However,  that  report  did  not  allow  for 
completely  arbitrary  pole  assignment,  and  the  modelling  assumptions  consider¬ 
ably  restricted  the  class  of  systems  for  which  it  was  applicable.  The 
results  given  in  (32)  overcome  these  problems. 

The  problem  of  finding  reduced  order  models  for  high  order  systems,  some¬ 
times  referred  to  as  the  "model  reduction  problem",  is  an  important  one  to  the 
practicing  engineer  since  it  is  difficult  to  apply  the  design  procedures  of 
modern  and  classical  control  theory  to  high  order  systems.  Nuoierous  solutions 
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have  been  proposed  during  the  past  two  decades  and  the  survey  of  Bosley  and 
Lees  (33)  represents  an  excellent  sumnary  of  the  available  methods.  A  number 
of  methods  are  based  on  first  deriving  transfer  functions  or  state-space 
models  for  the  high  order  systems  and  then  simplifying  these  models,  e.g.  (34), 
(35),  (36).  Other  methods  use  time  or  frequency  response  data  to  directly  fit 
low  order  models  (37),  (38),  (39).  Our  procedure  falls  into  this  latter  cate¬ 
gory,  since  it  determines  a  model  which  matches  the  frequency  response  of  the 
original  high  order  system  at  a  certain  set  of  prespecified  frequencies.  Its 
primary  advantage  lies  in  the  simplicity  of  implementation.  In  particular,  no 
intermediary  high  order  model  need  be  calculated,  only  one  test  input  need  be 
used,  and  the  calculation  of  model  parameters  only  requires  the  solution  of  a 
simple  set  of  linear  equations.  The  model  parameters  can  also  be  obtained  as 
the  output  of  an  analog  adaptive  network,  since  the  algorithm  makes  use  of  the 
generalized  equation  error  identification  scheme  proposed  by  Lion  (40).  Most 
importantly,  our  algorithm  readily  generalizes  to  the  multiple  input-multiple 
output  (multivariable)  case  where  classical  frequency  and  time  domain  proced¬ 
ures  become  cumbersome  to  apply. 

More  specifically,  in  (11),  we  present  a  hew  method  for  obtaining  simple 
low  order  models  whose  dynamical  behavior  approximates  that  of  more  complex, 
higher  order,  stable  linear  systems.  The  low  order  model  is  determined  by 
applying  an  identification  procedure  to  input-output  data  obtained  by  driving 
the  original  system  with  a  special  periodic  test  signal.  We  prove  that  in  the 
scalar  case  a  Lion-type  model  adjustiaent  identifier  will  determine  a  constant 
kfh  order  model  of  an  n^^  order  (k  <n)  system  provided  the  system  input  con¬ 
sists  of  exactly  k  distinct  sinusoids.  This  order  model  will  approximate 
the  higher  order  systems  in  the  sense  that  its  frequency  response  matches  that 
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of  the  model  at  the  k  input  frequencies  (provided  the  model  obtained  is 
stable).  We  then  show  that  this  result  can  be  extended  in  a  very  natural  way 
to  the  multivariable  case.  We  finally  demonstrate  by  example  that  this  pro¬ 
cedure  can  produce  excellent  low  order  models  when  such  models  exist. 

6.  Multi-Purpose  Controllers 

Frequency  domain  mothods  have  always  dominated  control  system  design  in 
the  scalar  (single  input/output)  case,  when  compared  to  the  more  "modern” 
state-space  or  differential  operator  methods,  due  to  the  relative  simplicity 
of  the  resulting  controllers  and  their  ability  to  function  acceptably  over  a 
rather  wide  range  of  plant  parameter  variations;  i.e.  their  robustness.  It  is 
not  surprising,  therefore,  that  numerous  studies  have  been  made  to  "extend" 
various  frequency  domain  techniques  to  the  multivariable  case  in  order  to 
simply  and  reliabley  achieve  a  diversity  of  desired  design  goals.  In  most 
cases,  however,  direct  extensions  of  scalar  frequency  dosiain  procedures,  such 
as  the  Nyquist  stability  criteria  or  the  root  locus,  are  not  possible  and 
often  rather  complex  modifications  have  to  be  made  to  existing  theories  in 
order  to  achieve  appropriate  design  objectives.  Further  complicating  the 
picture  is  the  fact  that  noninteraction  (or  decoupling)  is  often  an  additional 
design  objective  in  the  multivariable  case,  and  a  completely  decoupled,  stable 
system  cannot  always  be  achieved  by  the  relatively  simple  feedforward  control¬ 
lers  obtained  by  multivariable,  frequency  domain  methods. 

On  the  other  hand,  the  so  called  "modern"  methods  which  have  generally 
relied  on  exact  knowledge  of  the  plant,  are  continually  being  improved  upon 
and  extended  to  take  into  consideration  parameter  uncertainty  and/or  varia¬ 
tions;  i.e.  robustness  is  becoming  increasingly  important  in  designs  based  on 
state-space  or  differential  operator  methods.  Although  these  "modern"  methods 
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generally  imply  more  complex  controller  configurations,  than  those  associated 
with  frequency  domain  methods,  they  are  less  heuristic  to  iaiplement  and  can 
generally  achieve  more  than  is  possible  with  the  simpler  controllers  designed 
by  frequency  domain  methods.  Moreover,  with  the  ever  increasing  utilization 
of  computers  in  the  control  loop,  it  may  be  argued  that  controller  simplicity 
is  no  longer  as  important  as  it  once  was,  and  one  might  therefore  expect  to 
see  more  complex  controllers  being  used  in  future  applications. 

In  light  of  these  observations,  a  new  procedure  has  been  developed  for 
designing  controllers  which  simultaneously  achieve  a  variety  of  desired  design 
goals  in  deterministic,  unity  feedback,  linear  multivariable  system.  This 
procedure  involves  an  appropriate  combination  of  work  delineated  in  three 
earlier  reports,  the  first  dealing  with  the  pole  placesient  stabilization  of 
stabilizable  systems  (41),  the  second  with  skew  prime  polynomial  matrices  and 
their  importance  in  extending  the  internal  model  principle  to  the  nnjltivar- 
iable  case  (42),  and  the  third  with  output  regulation  and  tracking  with  zero 
steady-state  error  (43).  All  three  of  these  efforts  are  appropriately 
combined  in  (12)  to  yield  a  new  algorithm  for  the  systematic  design  of  a 
"three  part"  sniltivariable  controller  which  simultaneously  insures 

(a)  a  non  interactive  or  decoupled  closed  loop  design, 
well  as  closed  loop  stability, 

(b)  complete  and  arbitrary  closed  loop  pole  placement,  which  implies 

desired  (single  loop)  transient  performance  as  well  as  closed  loop 

stability, 

(c)  zero  steady-state  errors  between  the  plant  outputs  and  any 
nondecreasing  deterministic  inputs. 
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(d)  complete  steady-state  output  rejection  of  nondecreasing  determin¬ 
istic  disturbances,  and 

(e)  robustness  with  respect  to  stability,  disturbance  rejection,  and 
zero  error  tracking  for  rather  substantial  plant  parameter 
variations. 

Our  development  employs  the  more  "modern"  (Laplace  transfonsed)  differ¬ 
ential  operator  approach  (27)  for  controller  synthesis,  which  involves 
transfer  matrix  factorations  and  the  manipulation  of  polynomial  siatrices 
in  the  Laplace  operator  s. 

Reference  (43)  (Remark  6  in  particular)  indicates  how  one  might 
translate,  to  differential  operator  form,  the  results  in  (44),  (45)  and  (46) 
in  order  to  design  a  robust  compensator  which  yields  either  regulation  or 
tracking.  These  results  are  extended  in  (12)  to  include  sisniltaneous  regul¬ 
ation  and  tracking  using  a  single  robust  "internal  model",  M**),  of  both  the 
disturbance  and  the  reference.  We  also  show  how  the  remainder  of  our  feed¬ 
forward  compensator  can  be  chosen,  along  with  the  decoupled  open  loop  poles, 
to  achieve  an  overall  set  of  desired  closed  loop  decoupled  poles  under  unity 
feedback.  Most  importantly,  we  have  successfully  "combined"  the  results  given 
in  the  earlier  references  to  present,  for  the  first  time,  a  formal  proof  of 
our  ability  to  achieve  all  of  the  desired  design  goals  simultaneously. 

7.  Constrained  Control 

In  the  control  of  linear  multivariable  systems,  constraints  are  often 
imposed  on  both  the  complexity  of  the  controller  as  well  as  its  placement 
relative  to  the  system's  inputs  and  outputs.  An  interesting  question  relative 
to  this  observation  is  what  changes  occur  to  the  dynamical  relations  between 


any  input/output  pair  aa  the  reault  of  applying  constant  gain  output  feedback 


between  any  (i-th)  output ^  y^and  any  (j-th)  input,  uj.  This  question  is 
clearly  related  to  that  of  decentralized  control;  i.e.  of  determining  the 
conditions  under  which  one  or  more  "local"  output  feedbacks  can  be  applied  to 
insure  complete  state  controllability  of  the  system  through  a  selected  input 
or  set  of  inputs.  The  solution  to  the  rather  general  question  which  is  posed 
in  (47)  contains  the  elements  for  resolving  the  decentralized  control 
question,  as  well  as  other  "constrained"  control  questions,  in  a  new  and 
efficient  manner. 

The  particular  treatment  employed  in  (47)  assumes  knowledge  of  the  (pxm) 
rational  transfer  matrix,  T(s),  which  characterizes  the  dynamical  behavior  of 
the  system  under  investigation.  By  then  employing  a  relatively  right  prime 
factorization,  R(s)P(s)~^,  of  T(s)  with  P(s)  in  unique  (upper  right  triangu¬ 
lar)  Hermite  form,  significant  new  insight  is  obtained  relative  to  the  changes 
which  occur  in  all  of  the  pxm  elements  of  T(s)  when  any  output  y£  is  fed  back 
to  any  input,  uj  through  a  constant  gain  element  gji« 

To  be  more  specific,  let  us  consider  a  linear  aniltivariable  system  whose 
dynamical  behavior  is  specified  by  a  (pxm)  proper,  rational  transfer  matrix 
T(s)  ■  [t£j(s)  ].  Let  us  further  assume  that  n  denotes  the  order  of  any 
minimal  (state-space)  realization  of  T(s),  and  that  all  n  poles  of  T(s)  are 
distinct  (such  an  assumption  serves  to  simplify  this  presentation,  but  is  not 


r  (.8) 

really  necessary  in  general).  If  each  rational  element  t£j(s)  ■  p'  of 

T(s)  is  then  represented  in  prim,  factored  form,  one  can  imnediately  deter¬ 
mine  th  "controllability/observability  properties"  of  each  of  the  n  modes  or 
poles  by  inspection  of  T(s).  To  illustrate,  if  p  "  m  *3,  and 
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T(s)  « 


1  s*l  1 

s+1  s+4  s+1 

JL  0 

s+2  s+4 

1  s+2  s+1 

s+3  (s+3)  (s+5)  s+6 


then  it  is  clear  the  TCs)  has  6(  *n)  distinct  modes  at  s  ■  -1,  -2,  -3,  -4,  -5 
and  -6.  Furthermore,  since  the  mode  at  s  ■  -1  appears  only  in  the  first  row 
and  first  and  third  columns  of  T(s),  it  is  observable  only  at  the  first  out¬ 
put,  yx  and  controllable  only  by  the  first  and  third  inputs,  ux  and  U3, 
respectively.  Similar  observations  can  readily  be  made  regarding  the  other 
(five)  system  modes;  e.g.  the  mode  at  s  >■  -4  is  controllable  only  by  U2  and 
observable  only  at  yx  and  y2. 

One  of  the  fundamental  questions  in  constrained  control  is  the  following 
What  system  modes  become  controllable  via  u;^  if  y£  is  fed  back  to  uj 
arbitrary  gain  gj£,  as  depicted  in  the  figure  below?  It  might  be  noted  that 
only  the  outputs  and  inputs  pertinent  to  our  discussion  are  depicted. 


With  respect  to  this  fundamental  decentralised  control  question,  we  have 


been  able  to  prove  the  following  important  new  result: 


Theorem:  If  t£]((s)  ■ 


u^(.) 


f  0,  and  yi(s)  is  fed  back  to  uj^|((s) 


via  gji,  then  all  modes  controllable  via  u;((s)  before  feedback  will  remain 
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controllable  via  ui((3)  irrespective  of  gjj  (although  certain  of  these 
controllable  modes  could  vary  with  gji).  In  addition,  all  modes  controllable 
via  uj,  but  not  controllable  via  uj^,  will  become  controllable  via  uj^  for 
almost  any  gj£  except  those  which  are  both  (i)  unobservable  at  yj  and  (ii) 


zeros  of  riit(s),  the  prime  transfer  function. 

We  now  note  that  in  view  of  Theorem  1  and  the  given  (3  x  3)  T(s)>  if  y3Bi. 
is  fed  back  to  u£i>j  via  gi3>  then  the  U23i(  controllable  modes  at  s  »  -4  and  s 
■  -5  will  remain  controllable  via  U2  as  will  the  "new"  U2  controllable  mode  at 
8  “  ~3+gi3.  In  addition,  the  uj  controllable  mode  at  s  =  -1  will  become  con¬ 
trollable  via  U2  for  almost  any  gi3.  However,  the  u^  controllable  mode  at  s  = 


-2  will  not  become  controllable  via  U2  after  feedback  because  it  is  (i)  unob- 

.  y-,(s)  .  - 

servable  at  y3  and  (ii)  a  zero  of  night  be 

noted  that  if  we  now  set  k  *  3,  our  theorem  clearly  indicates  that  the  U3  ^0^. 


trollable  modes  at  s  =-l  and  -6  will  remain  U3  controllable.  In  addition, 
the  ui  controllable  modes  at  -2,  and  -3  +  gi3  will  become  controllable  via  U3 
for  almost  any  g£3  after  to  uj  feedback. 

We  merely  remark  here  that  a  "dual"  result  involving  yj^  observable  modes 
after  y£  to  uj  feedback  via  gj£  has  also  been  formally  stated  and  established. 
Furthermore,  when  combined,  these  two  results  provide  a  direct  algorithm  for 
determining  the  changes  which  occur  in  all  (pxm)  elements  of  the  closed  loop 
transfer  matrix  of  the  system.  Such  results  can  also  be  used  to  develop  new 
control  techniques  for  multivariable  systems  under  decentralized  control  con¬ 
straints.  We  are  currently  working  on  the  development  of  such  techniques. (48) 
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It  has  long  been  of  interest  to  find  explicit  conditions  for  a  real, 

canonical  linear  system  (Cpxnt Anxnt^nxm^  to  be  completely  assignable;  i.e.  to 

have  the  property  that  for  each  real,  monic,  nth  degree  polynomial 

n  j 

a(X)  =  x”  +  Z  a  X  (A) 

i=l 

there  exists  a  real,  constant  output-feedback  matrix  F^ixp  for  which 

det(Xl-A-BFC)  =  a(  X)  (B) 

It  has  also  been  of  interest  to  identify  those  systems  which  are  'generically 
assignable':  (C,A,B)  is  generically  assignable  if  the  set  of  all  coefficient 
vectors  (ax ,a2, • • . ,a{^) '  associated  with  (A)  for  which  thee  exist  real  F 
satisfying  (B)  is  open  and  dense  in  R". 

Apart  from  the  cases  when  either  C  has  independent  columns  or  B  has  inde¬ 
pendent  rows  (which  can  be  dealt  with  using  state-feedback  theory),  little  is 
known  about  either  complete  or  generic  assignability.  Perhaps  the  sharpest 
result  to  date,  due  to  Kimura  (49)  and  others,  asserts  that  for 
'almost  every'  linear  system,  generic  pole  assignment  is  possible  provided 
n  ^  m  +  p  -  1.  Willems  and  Hesselink  (50)  show  that  for  almost  every  system 
with  m  =  p  =>  2  and  n  *  4,  generic  pole  assignment  is  not  possible,  even 

though  the  number  of  free  parameters  in  F  equals  n.  Using  a  vesion  of  the 
(51)  implicit  function  theorem,  Hermann  and  Martin  prove  that  for  almost  every 
linear  system  whose  first  n  Markov  matrices,  CB,CAB, . . . ,^~^B  are  linearly 
independent,  generic  pole  assignment  is  possible  provided  F  is  allowed  to  be 
complex-valued.  Brockett  and  Byrnes  (52)  take  advantage  of  certain  classical 
ideas  based  on  elimination  theory  to  devleop  a  formula  which  gives  values  of 
m,n  and  p,  for  which  almost  every  linear  system  is  generically  assignable  -but 
their  developpient  is  not  constructive. 
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In  (53),  we  develop  a  useful  foraulacion  of  the  assignment  problem  for  p 
»  2  and  provide  constructive  solutions  for  some  special  cases.  For  m  *  2  and 
n  ■  A  we  characterize  the  classes  of  genetically  assignable  and  completely 
assignable  systems.  We  show  that  for  m  *  3  and  n  **  6,  almost  every  linear 
system  is  genetically  assignable;  we  do  this  by  reducing  the  assignment 
problem  to  the  problem  of  finding  a  real  root  of  a  real  5th  degree  polynomial 
in  one  variable.  That  such  a  polynomial  should  exist  is  in  agreement  with 
the  Brockatt-Byrnes  formula  derived  in  (52). 

8.  Digital  Control  of  Multivariable  Systems 

Perhaps  the  most  important  new  issue  in  multivariable  control,  is  the 
question  of  actual  implementation  of  the  final  design.  Whereas  ten  years  ago 
many  of  the  designs  were  too  complex  and  expensive  to  implement,  microproces¬ 
sor  technology  offers  an  exceptional  opportunity  for  overcoming  this  limita¬ 
tion.  Currently,  there  is  a  revolution  in  what  sorts  of  controller  designs 
can  and  should  be  physically  employed  in  practice.  Because  of  the  advent  of 
low  cost  microprocessors,  future  designs  will  be  not  only  more  reliable,  but 
also  less  expensive.  The  proper  and  effective  utilization  of  microprocessors 
for  the  control  of  complex  multivariable  processes  is  a  new  and  exciting 
research  area  which  presents  many  challenges  to  the  control  engineer. 

It  may  be  noted  at  this  time  that  much  as  been  done  and  is  known  regard¬ 
ing  discrete  control,  expecially  in  the  scalar  case,  i.e.,  texts  and  formal 
articles  dealing  with  sampled  data  systems  have  been  commonplace  for  years  and 


are  too  numerous  to  mention.  In  most  of  these  texts  and  articles  however, 


microprocessors  are  not  explicitly  mentioned  as  the  discrete  time  device  and, 
perhaps  more -importantly,  the  limited  memory,  wordlength,  and  speed  associated 
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with  microprocessors  are  not  usually  taken  into  consideration.  Most  of  the 
current  theory  as  well  as  the  applications  of  digital  control  still  assume  a 
rather  simple  single  input/output  system  to  be  controlled,  usually  by  a  dis¬ 
crete  device  with  sufficient  memory,  wordlength,  and  speed  to  approximate  that 
which  would  be  attainable  using  a  continuous  controller.  Nonetheless,  micro¬ 
processor  technology  has  given  renewed  impetus  to  the  area  of  discrete 
control,  and  contemporary  texts  and  articles  in  this  field  are  becoming 
increasingly  aware  of  the  potential  as  well  as  the  limitations  of 
microprocessors . 

With  respect  to  microprocessor  control  of  complex  multivariable  systems, 
much  theoretical  as  well  as  practical  work  remains  to  be  done.  The  work 
already  accomplished  in  this  area  involves  the  first  complete  solution  to 
deadbeat  error  tracking  in  discrete  multivariable  systems  (54).  In  particular, 
by  employing  a  polynomial  matrix  representation  for  the  dynamical  behavior  of 
a  discrete  system  in  the  delay  operator  d  (which  is  ^-transformable  to  z~^), 
the  necessary  and  sufficient  condition  has  been  found  for  achieving  a  discrete 
error  signal  between  reference  input  and  plant  output.  The  condition  involves 
the  skew  primeness  of  two  polynomial  matrices  in  d,  and  extends  the  internal 
model  principle  to  the  discrete  case.  Work  is  proceeding  in  actually  applying 
the  theory  to  the  microprocessor  control  of  a  linear  motor  in  our  new 
Laboratory  for  Engineering  Man/Machine  Systems  (LEMS). 
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